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ABSTRACT 

We present intermediate resolution, wide wavelength coverage spectra for a complete 
sample of 14 compact radio sources taken with the aim of investigating the impact of 
the nuclear activity on the circumnuclear interstellar medium (ISM) in the early stages 
of radio source evolution. We observe spatially extended line emission (up to ~20 kpc) 
in the majority of sources which is consistent with a quiescent halo. In the nuclear 
apertures we observe broad, highly complex emission line profiles. Multiple Gaussian 
modelling of the [O IIIJA5007 line reveals between 2 and 4 components which can have 
velocity widths (FWHM) and blueshifts relative to the halo of up to ^2000km s^^. We 
interpret these broad, blueshifted components as material in outflow and discuss the 
kinematical evidence for jet-driven outflows as previously proposed for PKS 1549-79 
and PKS 1345-1-12. Comparisons with samples in the literature show that compact 
radio sources harbour more extreme nuclear kinematics than their extended counter- 
parts, a trend seen within our sample with larger velocities in the smaller sources. 
The observed velocities are also likely to be influenced by source orientation with 
respect to the observer's line of sight. Nine sources have associated HI absorption. 
In common with the optical emission line gas, the HI profiles are often highly com- 
plex with the majority of the detected components significantly blueshifted, tracing 
outflows in the neutral gas. The sample has been tested for str atiflcation i n the ISM 
(FWHM/ionisation potential/critical density) as suggested by Holt ct al. (2003) for 
PKS 1345-1-12 but we find no significant trends within the sample using a Spearman 
Rank analysis. This study supports the idea that compact radio sources are young 
radio loud AGN observed during the early stages of their evolution and currently 
shedding their natal cocoons through extreme circumnuclear outfiows. 

Key words: ISM: jets and outfiows - ISM: kinematics and dynamics - galaxies: active 
- galaxies: ISM - galaxies: kinematics and dynamics 



1 INTRODUCTION 

We currently know little about the early evolution of power- 
ful extragalactic radio sources, in particular, how the pres- 
ence of a radio-loud AGN influences the evolution of the host 
galaxy, and how the ISM affects the expansion of the radio 
jets. Recent developments in both observation and theory 
have shown the importance of AGN feedback in galaxy evo- 
lution. For example, there exist close correlations between 
the mass of the central black hole and the properties of the 
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galaxy bulge (e.g. Ferrarese fc Merrittll2000l; Gebhardt et al.l 



I2OO0I : iTremaine et all |2002| : iMarconi fc Huntll2003l). T heo- 
retical analys es by, for example. ISilk fc Reed (|l998l ) and 
iFabianI (Il999h. and more rece ntly, simulations by. for ex - 
ample, iDi Matteo et all (|2005l ) and iHopkins et~al] l|2005l ). 
describe how, after reloc ating to the centre o f the galaxy 
remnant after a merger l|Heckman et al.lll986l ). the super- 
massive black hole grows through accretion, becomes active 
(a proto-quasar) and then sheds its enshrouding cocoon (de- 
posited by the merger) through outflow s driven by powerful 
quasar winds l|Balsara fc Krolick|[l993l ') . With time, the cen- 
tral regions will be cleared of fuel, starving both the central 
AGN and any star forming regions in the the bulge and 
activity will cease. 
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Clearly, understanding AGN-induced feedback is of vi- 
tal importance for understanding the evolution of both AGN 
and their host galaxies. However, there remain many un- 
certainties in the models due to a lack of observational evi- 
dence, and feedback is often inserted as a 'black box'. Whilst 
the models described above provide a good theoretical de- 
scription of the shedding of the natal cocoon during the 
early evolutionary stages, these scenarios often assume that 
the feedback is dominated by quasar- induced winds. Whilst 
this may be true for radio-quiet AGN, in radio-loud AGN, 
the expanding radio jets may also provide a significant con- 
tribution to the over all feedback of the AGN through jet- 
induced outflows (e.g. lBicknell fc Dopita|[l997l ). particularly 
during the early stages of the radio source evolution when 
the radio jets are still on the same scales as the natal co- 
coon. Such a contribution on small scales should not be sur- 
prising given the increasing support for significant feedback 
from extended radio sources, for example , in the quench- 
ing of cooling flow s in massive haloes (e.g. iBest et ahlbOOfil : 
ICroton et aiTl2006l : iBwyer et al.ll2006l ). 

With this in mind, compact radio sources, compris- 
ing the Gigahertz-Peaked Spectrum radio sources (GPS: 
D < 1 kpc) and the larger Compact Steep Spectrum ra- 
dio sources (CSS: D < 15 kpc), form a large proportion 
of the bright (centimetre - wavelength-s elected) radio source 
population ('-^40%; e.g. iO'Deaj Il998t ) particularly suited 
to studying AGN-induced feedback. First, GPS and CSS 
sources are cu rrently believed t o be compact due to evolu- 
tionary stage (iFanti et al.lll995l ) with estimates of dynam- 
ical and radio spectral ages of tdvn ^ 1 0^ — 10^ yr (e.g. 
lOwsianflt et~a r 1998: Tschager et al.|[200ol ) and tsp < lO"* yr 
fe.g. iMurgiret al. 19991 7 respectively. This is in preference 
to the frustration scenario in which the ISM is so dense that 
the radio jets cannot escape, and the radio so urce remain s 
confi ned and frustrated for its entire lifetime ( van Breug el 
Il984l ). Second, the small-scale radio jets will be on the same 
scales as the natal cocoon of gas and dust and so will readily 
interact with it. Finally, radio- loud sources contain all po- 
tential outflow driving mechanisms (quasar-induced winds, 
radio-jets and starburst-driven superwinds) and so are the 
only objects in which the relative importance of the dif- 
ferent effects can be assessed in individual objects (e.g. 
iBatcheldor et al.ll2007^ . 

The most direct way to probe the kinematics and phys- 
ical conditions of the nuclear gas is through high quality 
optical spectroscopy. However, to date, much of the work on 
compact radio sources has either concentrated on the radio 
wavelength region or has relied on low signal to noise optical 
spectra. Despite this, an early optical spectroscopic study of 
a large sample of compact radio sources suggested evidence 
for non-gravitational motions (indicative of flows of gas) in 
the form of highly broadened emission line profllcs which are 
often asymmetric with blue wings (Gcldcrman & Whittle 
1 19941 ). However, the potential difi'erences between the spec- 
tra of compact and extended sources were not quantified in 
this early study. 

With new, deeper spectra, it becomes clear that the un- 
usually broad nuclear emission lines are due to highly com- 
plex emission line profiles requiring multiple Gaussians to 
model them. The first concrete evidence for fast outflows in 
the optical e mission line gas in a co mpact radio source was 
reported by iTadhunter et al.] (|200lh in the southern com- 



pact flat spectrum radio source PKS 1549-79. From their 
low resolution optical spectra, the high ionisation emission 
lines (e.g. [O HI]) were both broader (FWHM ~ 1350 km 
s~^ compared to ~ 650 km s~^) and hlueshifted by ~ 600 
km s~^ with respect to the low ionisation lines (e.g. [O II]). 
Recent higher resolution follow up on this source has shown 
that the outflowing compo nent is present in all lines with a 
velocity o f 679 ± 20 kms~^ l|Holt et al.l2006l : hereafter H06). 
Similarly. iHolt et 31.1(120031 ) (hereafter H03) reported a much 
more extreme outflow (up to 2000km s~^) in the GPS source 
PKS 1345-1-12. In both sources, the outflows are believed to 
be driven by the small scale radio jets expanding through 
dens e circumnuclear cocoon s of gas and dust (see e.g. Figure 
2 in ITadhunter et al.ll200ll and Figure 4 in H03). A recent 
high-resolution optical and radio imaging study (HST and 
VLBI) was able to confidently rule out large-scale starbursts 
driven winds, but failed to distin guish between jet-drive n 
outfiows and AGN-induced winds l|Batcheldor et aLlbOOTj ). 
However, previous HST imaging studies of compact radio 
sources have suggested evidence for an alignment effect on 
small scales in some CSS sources l|de Vries et al.lll997l . [l999l : 
Ax on et al. ,2000. ) which would support the jet-driven out- 
flows scenario. For PKS 1345-1-12, due to the differences in 
velocity width between the emission line components in dif- 
ferent lines, H03 suggested a stratifled ISM which may be 
the result of gradients in density and/or ionisation potential 
across the circumnuclear ISM. 

Hence, we have obtained intermediate resolution (4- 
6^ optical spectra with good signal-to-noise over a large 
spectral range, with the WHT, NTT and VLT, in order to 
search for such outflows in a complete sample of 14 compact 
radio sources. In this paper, the first of two, we present the 
kinematical data for the whole sample in order to search for 
the signatures of outfiows in the nuclear emission line gas. 
In sections 2 and 3, we describe the sample selection, ob- 
servations and data reduction procedure. In section 4, we 
model the emission lines, both in the extended halo (to de- 
termine the rest frame of each system) and in the nucleus 
(to search for outfiows). Our data are of sufficient resolu- 
tion to model the lines with multiple Gaussian components. 
In addition, in section 5 we have modelled the lines with 
single Gaussians to allow comparisons between the sample 
presented here and other samples (compact and extended) 
in the literature, and to investigate the possibility of density 
and/or ionisation gradients in the ISM in some sources. In 
this paper we concentrate on the emission line kinematics 
and the ionisation of the emission line gas will be discussed 
in a second paper. 

Throughout this paper, we assume the following cos- 
mology: Ho = 71 km s"\ fio = 0.27 and = 0.73. 



2 SAMPLE SELECTION 

Our complete sample comprises 14 compact radio 
sources split across the northern and southern hemispheres. 
In the north, we have observed a complete sample of five 
sources (3 CSS, 2 GPS wfth z < 0.4, 12*^ < RA < IS'' 

^ This corresponds to r^200-300km s"'^, dependent on redshift, 
for [O III]5007A for the objects in our sample. 
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Table 1. Properties of the sample: (a) radio source; (b) optical ID where G = galaxy, Q = QSO; (c) radio source type where CSS = 
Compact Steep Spectrum, GPS = Gigahertz-Peaked Spectrum, CF = Compact Flat spectrum and CC = Compact Core; (d) redshift * 
redshifts from the literature, for the accurate redshifts derived in this study, se e Table \S\ (e) radio spe ctral index where F^, oc and 
for ol'^gfj^ (for all but 3C 213.1, 3C 277.1 and 3C 303.1, these are taken fro m lWall fc Peacockl lIlQSSi '); for the remaining three sources 
a is calculated from the flux densities presented in iKellermann et al.l lll969l) ); (/) radio luminosity at 5GHz; (g) HI 21cm absorption 
velocity (km s~^); (h) HI 21cm absorption FWHM (km s~^); (i) Angular size of radio source in milli arcseconds; (j) Projected linear 
size of radio source in parsecs; (fc) Radio axis position angle. For references, see the detailed biographies of each sources in the text, f 
tentative HI detection with optical depth, t = 0.003 (Morganti 2004, private communication) where the central velocity is assumed to 
be the midpoint of the entire absorption profile. " taken from the literature and may be inconsistent with the sizes in column (j), ^ 
compact core co mponent in 3C 213.1 a nd 3C 459, extended radi o sourc e component in 3C 213.1 and 3C 459. For 3C 459 the different 
HI results are: ™ iMorganti et al.l l|20o3), ^ IVermeulen et al.l ||2003|) and s ICupta et all (.2006, '). 



Object 


ID 


Type" 




a 


Radio 


HI 21cm absorption 


Angular 


Linear 


Radio 












luminosity 


^ helio 


FWHM 


size 


size 


PA 










1 


log P5GH. (W Hz-1) 


(km s"-*-) 


(km s~-^) 


(arcsec) 


(kpc) 


(°) 


(a) 


(b) 


(c) 


id) 


(e) 


(/) 


(ff) 


(h) 


(i) 


(J) 


(k) 


3C 213.1 


G 


CC 


0.195 


0.30 


26.41 


0.19395 


115 


6.0* 


19.1* 


-61? 


















40= 


128^= 


-20 


3C 268.3 


G 


CSS 


0.371 


0.88 


27.04 


0.37186 


101 


1.36 


6.9 


-15 














0.37227 


19 








3C 277.1 


Q 


CSS 


0.320 


0.61 


26.86 






1.67 


7.7 


-49 


4C 32.44 


G 


GPS 


0.369 


0.60 


27.30 


0.36843 


229 


0.06 


0.31 


-50 


PKS 1345+12 


G 


GPS 


0.122 


0.44 


26.69 


0.12184 


150 


0.15 


0.33 


160 


3C 303.1 


G 


CSS 


0.267 


1.08 


26.43 






1.80 


7.3 


-47 


PKS 0023-26 


G 


CSS 


0.322 


0.7 


27.43 


0.31890 


126 


0.68 


3.2 


-34 














0.32142 


39 








PKS 0252-71 


G 


CSS 


0.566 


1.14 


27.55 






0.24 


1.6 


7 


PKS 1306-09 


G 


CSS 


0.464 


0.65 


27.39 


0.46492t 


350t 


0.46 


15.7 


-41 


PKS 1549-79 


G 


CF 


0.152 


0.18 


27.00 


0.15224 




0.12 


0.3 


90 


PKS 1814-63 


G 


CSS" 


0.063 


0.91 


26.54 


0.06450 




0.41 


0.5 


-20 


PKS 1934-63 


G 


GPS" 


0.183 


0.88 


27.31 


0.18282 


18 


0.70 


2.1 


89 


PKS 2135-20 


G 


CSS 


0.635 


0.82 


27.58 






0.25 


1.7 


52 


PKS 2314-1-03 


G 


CC 


0.220 


0.97 


27.65 


0.21850™ 


-400 


0.20* 


0.7* 


95 














0.21914^ 


130 


10.00^= 


35.2^= 





0.21837S 71 
0.21839S 121 
0.21963S 164 



and 6 > 15°) from the subset of compact radio sources 
in the 3rd and 4th Cambridge Radio Catalogues (3C and 
4C respectiv ely, at 178 MHz) listed in Table 1 in the 
lO'Deal 11998') review, which is a master list of the CSS 
sample from Fant i et al. ( 19901) and the GPS sample from 
IStanghellini et aU (Il990i . 1 19961 ). A further seven sources (5 
CSS, 1 GPS, 1 compact flat-spectrum) form a complete red- 
shift (z < 0.7) and RA (12'' < RA < 5'') l imited subset of 
the 2 Jy sample of iTadhunter et "ahl (|l993l), which in turn 
is a su bset {z < 0.7, S < +10°) of the 2Jv lWaU fc Peacockl 
l|l985l ) radio survey at 2.7 GHz. We have also observed a 
further two radio sources (3C 213.1 and 3C 459) which are 
classified as compact core radio sources. Such sources have 
both extended radio structure (jets and lobes which are of- 
ten disconnected from the radio core) and a second, compact 
radio structure (complete with lobes), which may be mis- 
aligned from the large-scale source. Whilst the larger-scale 
structures are consistent with old, evolved radio sources, the 
inner, compact structures may represent a younger, more 
recently triggered active phase, and so these sources can 
be studied in conjunction with a sample of compact radio 
sources. 



The whole sample has a radio power range of 26 < 
log PsGHz (W Hz~^) < 28. Details of the sample are pre- 
sented in Table [2] Hence, the total sample includes 8 CSS, 
3 GPS, 1 compact flat-spectrum source and 2 compact core 
radio sources. The sample also includes both galaxy (13) 
and quasar (1) host morphologies. 

The subsamples were chosen to fulflU the following 
criteria: i) redshifts low enough to include [O III]A5007 in 
the spectra and similar for the subsamples; ii) RA and Dec 
ranges to be observed during the allocated runs and; iii) 
the samples were derived from well-studied samples (2Jy, 
3C & 4C) to enable good comparisons between compact 
and extended sources. 



3 OBSERVATIONS, DATA REDUCTION AND 
ANALYSIS TECHNIQUES 

Long-slit optical spectroscopic observations of the full sam- 
ple were obtained during two observing runs. The north- 
ern sample was observed on 12-14 May 2001 with ISIS, 
the dual arm spectrograph on the 4.2m William Herschel 
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Table 2. Log of observations. Column (c) denotes the arm/wavelength range of the spectrograph, (h) is the instrumental seeing 
(t indicates seeing measured using the object continuum - for all other sources, a star on the chip was used). Note, the values 
quoted in column (h) represent upper limits on the true seeing as 1) where the seeing was measured using the target galaxy, these 
measurements may be affected by the extended galaxy haloes i.e. they are not stellar; and 2) the measurements based on the profiles 

of stars along the slit only provide a good seeing estimate if the slit passes through the centres of the stellar image, otherwise the 
FWHM (and therefore seeing) will be over-estimated, (i) gives the width of the extracted nuclear aperture in arcseconds, centred 
on the eentroid of the nuclear continuum emission (quoted once per object in the table but the same aperture was used for both 
red and blue spectral ranges) and (j) denotes the photometric conditions (p = photometric, v = variable transparency). ★ denotes 
spectra taken aligned within 10 degrees of the radio source PA. 



Date 


Object 


'Arm' 


Exposure 




Slit 


A range 


Seeing 


Nuclear 


Notes 










PA 


width 






aperture 










(s) 


(°) 


(arcsec) 


(A) 


(arcsec) 


(arcsec) 




(a) 


(6) 


(c) 




(e) 


{/) 


(9) 


(h) 


(i) 


U) 



WHT/ISIS observations 



14/05/2001 


3C 213.1 


R 


3*1200 


80* 


1.3 


7700-9200 


1.5 ± O.lf 


2.5 


V 






B 


3^^1200 


80* 


1.3 


3300-6800 


1.8 ± O.lf 




V 


13/05/2001 


3C 268.3 


R 


2*1200 


155* 


1.3 


7850-9350 


1.4 ± O.Sf 


2.2 


p 






R 


3'^1200 


155* 


1.3 


6200-7700 


1.9 ± 0.2t 




p 






B 


5*1200 


155* 


1.3 


3700-6800 


2.3 ± 0.2| 




p 


13/05/2001 


3C 277.1 


R 


1*1200 


129* 


1.3 


7850-9350 


0.9 ± O.lf 


1.8 


p 






R 


2*1200 


129* 


1.3 


6200-7700 


0.9 ± O.lf 




p 






B 


3*1200 


129* 


1.3 


3700-6800 


1.4 ± O.lf 




p 


14/05/2001 


4C 32.44 


R 


3 1200 


105 


1.3 


7700-9200 


1.5 ± 0.2t 


2.5 


V 






R 


3*1200 


105 


1.3 


6200-7700 


1.4 ± 0.2t 




V 


14/05/2001 




B 


3*1200 


105 


1.3 


3300-6800 


1.9 ± 0.2t 




V 


12/05/2001 


PKS 1345+12 


R 


1 900 


104 


1.3 


62U0-7700 


1.3 ± 0.2 


2.2 


p 




(4C 12.50) 


B 


1*900 


104 


1.3 


3300-6800 


1.3 ± 0.2 




p 






R 


3*1200 


160* 


1.3 


6200-7700 


1.3 ± 0.2 




p 






B 


3*1200 


160* 


1.3 


3300-6800 


1.3 ± 0.2 




p 


14/05/2001 


PKS 1345+12 


R 


3*1200 


230 


1.3 


6200-7700 


1.7 ± 0.2 


2.2 


V 






B 


3*1200 


230 


1.3 


3300-6800 


1.7 ± 0.2 




V 


12/05/2001 


3C 303.1 


R 


5*1200 


130* 


1.3 


7700-9200 


1.8 ± O.lf 


2.5 


p 






B 


5*1200 


130* 


1.3 


3300-6800 


1.8 ± O.lf 




p 


NTT/EMMI observations 


















13/07/2002 


PKS 0023-26 


B 


3*1200 


-105 


1.5 


3700-7050 


2.5 ± 0.2t 


2.0 


p 






R 


2*1200 


-105 


1.5 


4400-11400 


2.3 ± 0.2t 




p 


13/07/2002 


PKS 0252-71 


R 


2*1200 


135* 


1.5 


5700-8700 


1.6 ± 0.1 


1.7 


p 






B 


2*1200 


135* 


1.5 


3700-7050 


1.5 ± 0.1 




p 


13/07/2002 


PKS 1306-09 


R 


3*1200 


135 


1.0 


5700-8700 


1.7 ± 0.2t 


2.3 


p 






B 


3*1200 


135 


1.0 


3700-7050 


1.8 ± 0.2t 




p 


12/07/2002 


PKS 1549-79 


R 


3*1200 


-5 


1.0 


5700-8700 


1.9 ± 0.1 




p 






B 


3*1200 


-5 


1.0 


3700-7050 


2.2 ± 0.1 




p 






R 


2*1200 


25 


1.5 


5700-8700 


0.8 ± 0.1 




p 






B 


2*1200 


25 


1.5 


3700-7050 


0.8 ± 0.1 




p 


12/07/2002 


PKS 1814-63 


R 


1*1200 


-72 


1.0 


5700-8700 


1.1 ± 0.2t 


1.0 


p 






R 


1*1200 


65 


1.0 


5700-8700 


1.3 ± 0.1 




p 






B 


1*1200 


65 


1.0 


3700-7050 


1.2 ± 0.1 




p 


12/07/2002 


PKS 1934-63 


R 


3*1200 


-20 


1.5 


5700-8700 


1.2 ± 0.1 


2.0 


p 






B 


3*1200 


-20 


1.5 


3700-7050 


1.2 ± 0.1 




p 


13/07/2002 


PKS 2135-20 


R 


3*1200 


-115 


1.5 


5700-8700 


1.4 ± 0.2t 


1.7 


p 






B 


3*1200 


-115 


1.5 


3700-7050 


1.4 ± 0.2t 




p 


12/07/2002 


3C 459 


R 


3*1200 


-175 


1.5 


5700-8700 


1.5 ± 0.2t 


1.7 


p 




(PKS 2314+03) 


B 


3*1200 


-175 


1.5 


3700-7050 


1.5 ± 0.2t 




p 






R 


1*1200 


95* 


1.5 


5700-8700 


1.5 ± 0.2t 




p 


VLT/FORS2 observations 


















24/09/2003 


PKS 1549-79 


R 


3*1200 


75 


1.3 


4950-8250 


2.0 ± 0.1 


1.5 


p 


24/09/2003 




B 


3*600 


75 


1.3 


3050-6000 


2.1 ± 0.1 




p 
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Telescope on La Palma. In the red, the TEK4 CCD was 
used with the R316R grating with 1x1 binning and two cen- 
tral wavelengths (~6950A and ~8450A, dependent on ob- 
ject). The resulting wavelength calibration accuracies, cal- 
culated using the standard error on the mean deviation of 
the night sky emission line wavelengths from published val- 
ues (o^erbrockelal] [l99i, were 0.06-O.llA (~6950A) and 
O.IO-O.53A (~8450A). The spectral resolutions, calculated 
using the widths of the night sky emission lines, were 3.3- 
3.7 ± O.lA (~6950A) and 3.8-4.1 ± O.lA (~8450A). In the 
blue, the EEV12 CCD was used with the R300B grating 
with 2x2 binning. The wavelength calibration accuracy was 
O.I-O.2A with a spectral resolution of 4.3-4.8 ± O.2A. The 
spatial scale was 0.36 arcsec/pixel. Further details, including 
the instrumental setups and the useful wavelength ranges of 
the spectra are summarised in Table 2. 

The southern sample was observed on 12-13 July 2002 
using the EMMI spectrograph on the ESO New Technol- 
ogy Telescope (NTT) on La Silla, Chile in RILD mode. The 
MIT/LL CCD was used with grisms #4, #5 and #6 to 
obtain spectra with central wavelengths 7935A, 5385A and 
7223A with 2x2 binning. The wavelength calibration accu- 
racies were 0.24A (7935A), 0.06-0. 15A (5385A) and 0.10- 
O.I7A (7223A) with spectral resolutions of 14 ± lA (7935A), 
5.6 ± O.lA (5385A) and 4.3-6.7 ± O.lA (7223A). The spa- 
tial scale was 0.33 arcsec/pixel. Again, further details are 
presented in Table 2. 

PKS 1549-79 was also observed with the F0RS2 spec- 
trograph on the ESO Very Large Telescope (VLT) on Cerro 
Paranal, Chile, in September 2003 to improve the S/N of 
the fainter features and resolve a seeing-slit width matching 
problem in the NTT data. The detailed analysis of the data 
for PKS 1345-hl2 and PKS 1549-79 have already bee n pre- 
sented in H03, H06 and [Rodriguez Zaurin et~aLl l|2007h . The 
spatial scale for the VLT observations was 0.25 arcsec/pixel. 

The aim of this project is to search for outflows 
in the circumnuclear regions using optical spectroscopy 
with sufficient spectral resolution and signal to noise to 
accurately separate and model the different components of 
the highly broadened and complex emission lines. In order 
to include all of the outflowing regions in the slit, and to 
ensure the spectra were of sufficient resolution, all spectra 
were taken with a 1.0-1.5 arcsec slit. To reduce the effects of 
differential atmospheric refraction, all exposures were taken 
at low airmass (sec z < 1.1) and/or with the slit aligned 
along the parallactic angle. Due to various observational 
constraints, we have only aligned the slit along the radio 
axis for approximately half of the sources, and the PAs are 
listed in Table 2. However, as this study is concerned with 
compact rather than extended radio sources, the feedback 
effects we are interested in will be confined to the nuclear 
regions and the mis-alignment of the slit with the radio axis 
will not affect our results. 



3.1 Data reduction and analysis 

The data were reduced in the usual way (bias subtraction, 
flat fielding, cosmic ray removal, wavelength calibration, flux 
calibration) using the standard packages in iraf. The two- 
dimensional spectra were also corrected for spatial distor- 
tions of the CCD. To reduce wavelength calibration errors 



due to flexure of the telescope and instrument (ISIS observa- 
tions), separate arcs were taken at each position on the sky. 
Such calibrations were not made for EMMI or FORS2 as the 
reported flexure is small. The calculated wavelength calibra- 
tion accuracies and spectral resolutions are given above. 

Comparison of several spectrophotometric standard 
stars taken with a wide slit (5 arcsec) throughout each run 
gave relative flux calibrations accurate to ±5 per cent. This 
accuracy was confirmed by good matching in the flux be- 
tween the red and blue spectra. The only exception was 3C 
303.1 and whilst this will pose problems for line ratios, the 
kinematic information is unaffected. In order to correct the 
spectra for atmospheric absorption features (e.g. A and B 
bands at ~7600 and ~6800 A, respectively), standard stars 
were observed with a narrow slit, matched to the slit width 
used to observe the objects. For the ISIS observations, these 
were observed at random during the run and the normalised 
spectra had to be scaled to match the depth of the absorp- 
tion features. For the NTT and VLT observations, a sepa- 
rate calibration star was observed for each object, close in 
both position and time to the target observations and so no 
scaling was required. The spectra of all sources were cor- 
rected for Gal a ctic e xtinction using the E(B -V) values from 
ISchlegel et all l|l998l ) and the lSeatonI (| 19791 ) extinction law. 

The spectra were extracted and analysed using the 
STARLINK packages FIGARO and dipso. 

3.2 Continuum subtraction 

Prior to emission line modelling of the nuclear apertures, 
the continuum was modelled and subtracted for most of 
the sources. Initially, a n ebular continuum wa s generated 
and subtracted following iDickson et aD (|l995l ) taking full 
account of reddening following H03. The remaining contin- 
uum was then modelled usi ng a customised idl minimum y'^ 
fitting programme (see e.g. lRobinsonll200ll : iTadhunter et al.l 
I2OO5I : iHolt et al.ll2007l for details). The modelling program 
allows up to three separate continuum components - Old 
(OSP, here set to 12.5 Gyr with no reddening) and Young 
SteUar P opulations (Y SP, reddened with < E(B-V) < 
1.6 using Seatoiil 19791 w ith age 0.01-5.0 Gyr) taken from 
iBruzual fc CharlotI l| 1991 ) and an AGN/power- aw compo- 
nent. We define the best fitting model as that with the least 
number of components required to adequately model both 
the overall SED and discrete stellar absorption features (e.g. 
Ca H-l-K, Balmer series). 

Continuum modelling and subtraction was necessary 
for the detection and accurate modelling of the broader 
emission line components. For completeness, details of the 
continuum models used for subtraction prior to emission 
line modelling are given in Table O Note, the continua 
of several galaxies in this sample have been more re- 
cently accurately modelled and we r efer readers to the cor- 
responding papers: PKS 134 5-1-12 (|Tadhunter et al.l l2005l : 
iRodrieuez Zaurin et al.1 120071 ): PKS 154 9-79 (H06); PK S 
0023-26, PKS 1549-79 and PKS 2135-2 (|Holt et al.1 1 20071 '): 
3C 213.1 and 3C 459 (| Wills et al.ll2007l ). 

The continua of three sources were not modelled. 3C 
277.1 has a pure quasar continuum (there is no evidence for 
an underlying stellar continuum), but it could not be repro- 
duced by a simple power-law, and so we have not modelled 
or subtracted the continuum. However, when modelling the 
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Table 3. Summary of the continuum modelling results (see also Section 3.2) for the nuclear aperture of all galaxies in the sample. 
The columns are: (a) object, (6) best fitting (lowest Xjoduced) model — OSP (12.5 Gyr) plus YSP (Gyr) + power-law components, (c) 
percentage (flux) contribution of OSP, (d) percentage (flux) contribution of YSP, (e) E(B-V)value (if any) used to redden the YSP, 
(/) percentage (flux) contribution of the power law component, {g) spectral index of the power-law component where F;^ oc A" , ih) 
minimum reduced of the model for the whole SED. The nuclear continuum was not modelled for three galaxies: * 3C 277.1: pure QSO 
spectrum; o PKS 1814-63: continuum is strongly contaminated by a nearby Galactic star which drowns out any continuum information 
from the galaxy; f 30 303.1: mis-matching between the blue and red arms as the blue arm was observed without a dichroic (due to the 
position of [O III] AA4959, 5007). For these three galaxies, see the text for a discussion of alternative techniques. More detailed studies of 
the stellar pop ulations in some of the radio sources in this sample ar e presented else where and we ref er readers to the following papers: 
PKS 1345+12 llTadhunter et ahllioosl: iRodrigucz Zaurin ct al.l l2007h . PKS 1549-79 l|Holt et al.ll200d ). PKS 0023-26 and PKS 2135-20 

0% 



iHolt et al.l | |2007^ and 30 213.1 and 30 459 Wills ct al. ( 2007 
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emission lines, care was taken to remove the BLR contami- 
nation. The radio source PKS 1814-63 lies close to a bright 
foreground Galactic star and so the spectra are strongly con- 
taminated by hght from this star. Fortunately our obser- 
vations were made during good seeing condition^ and we 
clearly resolve PKS 1814-63 from the star. Whilst it was not 
possible to model the continuum, for accurate modelling of 
the broader emission line components, we have attempted to 
remove the stellar spectrum. Because the star is point-like, 
its spectrum will have no significant spatial variation and the 
star itself could be used for subtraction. Hence, each frame 
was copied, inverted and then shifted in the spatial direction, 
so the centroid of the continuum was aligned. The inverted 
spectrum was then subtracted from the original spectrum 
to remove the contamination from the continuum emission. 
Finally, no corrections were made for 3C 303.1 due to the 
mis-matching of the blue and red arm spectra. 

3.3 Kinematic component definitions 

The emission lines and their components vary significantly 
in width (FWHM). Hence, for our analysis, we define the fol- 
lowing kinematical components for the Narrow Line Region 
(NLR) as used in H03: 

• narrow: FWHM < 600 km s~^ 

^ Table 2 gives the seeing measured from a star on the 2- 
dimensional spectra whereas the seeing measured using several 
stars on the acquisition image was 0.8 arcsec. 



• intermediate: 600 < FWHM < 1400 km s"^ 

• broad: 1400 < FWHM < 2000 km s"^; 

• very broad; FWHM > 2000 km s~^ 



4 RESULTS 

In this section we discuss the kinematic results derived from 
two-dimensional optical spectra. Initially, the extended line 
emission is used to establish the systemic velocities followed 
by detailed modelling of the emission lines in the nuclear 
apertures. 

4.1 Kinematics of the extended gaseous halo 

In order to investigate the impact of the activity on the cir- 
cumnuclear material, it is first necessary to establish the ex- 
act redshift of the galaxy rest frame. Many of the objects in 
the sample show clear evidence for extended emission lines, 
most notably the [O II]AA3727 doublet and [O III]A5007, 
and in some cases Ha A6562 (see below) . By modelling these, 
pixel by pixel, across the entire spatial extent, one can deter- 
mine the radial velocity profile of the host galaxy (see also 
e.g. H03,H06). All such spatial fitting was done before mod- 
elling and subtracting the continuum and the emission lines 
were fitted using free fits. We define 'free-fitting' as when 
the only constraints used for the emission line modelling are 
those set by atomic physics (i.e. the separation, ratio of in- 
tensities and velocity widths of line components in a doublet 
emitted by a single ion). The radial velocity profiles for all 
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Figure 1. Rest frame velocity profiles for all sources in the sample. The numerous small points represent the pixel- by- 
pixel fitting where: triangles (A,V) are [O IIJAA3727 and squares/diamonds (□,(>) are [O III] AA4959, 5007. In general, 
open symbols (A,V,n,0) are narrow components (where two narrow components are present, the second is traced 
by inverted triangles and diamonds), filled symbols represent the intermediate (A, I) and broad (T,^) components. 
Exceptions are: 4C 32.44 (A, I: broad component; T,^: very broad component), PKS 1306-09 (open and filled points 
are narrow components), PKS 1814-63 ([O III] and Ha emission plotted where □,■: [O III] and O,*: Ho), PKS 1934- 
63 & 3C 459: (open: narrow, filled: broad), PKS 2135-20 (open: intermediate, filled: broad). Stars (where shown) 
mark the velocity of the detected HI absorption. Sources to note are: PKS 0023-26 (two components, small: narrow, 
large: broad), PKS 1306-09 (vertical dashed line marks the range of velocities over which HI absorption has been 
tentatively detected, Morganti, private communication), PKS 1814-63 (in addition to the deep HI absorption line 
(star) , the range of vel ocities covered by the b road shallow absorption featur e are marked by the ver tical dashed line) 
and 3C 459 (large star: lMorganti et al. 1*2001^ detection, large inverted star: IVermeulen et al.l |2003) detection, small 
stars: three components detected by Gupta et al. (2006)). The large symbols (n,^),^,^) represent the components 
of [O HI] AA4959,5007 detected in the nuclear aperture where the symbols (open/filled etc) are consistent with the 
smaller points in the same plot. The large symbols and their horizontal error bars show the position and size of the 
nuclear aperture. Finally, for 3C 459 (PA 95), two further points mark an extended region of line splitting ([O III]; 
O: narrow, •: broad) All extended apertures will be discussed in a future paper. The horizontal dashed line in each 
plot indicates the assumed systemic velocity (derived from [O III]). Note, for PKS 1345-1-12, this figure includes 
all three PAs observed and is taken directly from H03. For PKS 1549- 79, we show the velocity profile along PA -5 
only, again, taken directly from H06. For more details see H03, H06 and lRodriguez Zaurin et al.l l|2007h . Note, where 
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Figure 1. Velocity profiles continued. 



sources are shown in Figure [3Ve have determined the 
systemic velocities using a variety of techniques: 

(i) Sources with resolved extended emission. Where the 
line emission line is clearly resolved we have assumed the 
extended narrow line emission traces the quiescent halo gas 
as in PKS 1345-P12. For two sources (3C 277.1 and PKS 
1934-63), the extended gas has settled into what appears to 
be a smooth rotation curve whilst for the remaining sources, 
the gas is more similar to that in PKS 1345-1-12 - consistent 



^ Note, for some objects, the velocity fields are only marginally 
resolved (see Section 3 and Table 2 for details of the spatial scale 
and the seeing) and this is why in some of the plots in Figure 1, 
the velocities do not appear to vary significantly. 



with gravitational motions in an elliptical gala^xy (velocity 
amplitude < 350 km s~^: iTadhunter et al.lll989r ) but not 
yet completely settled. This technique is the most reliable 
for determining the systemic velocity and we have used this 
for most sources in the sample. 

(ii) Sources with two narrow components. In four sources 
(3C 268.3, 4C 32.44, PKS 0252-71 and PKS 1306-09) we 
detect two narrow components in the nuclear aperture (for 
3C 268.3, and marginally PKS 1306-09, these components 
are also extended and spatially resolved). Free fitting of the 
emission lines gives similar but slightly different FWIfM. 
Forcing the two narrow narrow components to have the 
same FWHM also gives a good fit, and so we assume the 
two narrow components are emitted by the same mecha- 
nism/structure. Hence, the narrow line splitting we observe 
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Figure 1. Velocity profiles continued. 



could represent a rotation curve which is unresolved, or trace 
a bi-polar outflow in the gas. For these sources, we can there- 
fore assume with reasonable confidence that the systemic 
velocity is located between the two narrow components. 

(iii) Other methods. For two sources, we have little confi- 
dence in determining the systemic velocity. For PKS 1814- 
63, the foreground star obscures our view of the galaxy. Us- 
ing the continuum subtraction technique described in Sec- 
tion 3.2, we have been able to model the extended line emis- 
sion along two PAs. Along PA 65 the shape of the profile 
is unclear whilst along PA -72 a more defined shape is ob- 
served, which may be consistent with half of the rotation 
curve. For PKS 2135-20, we are the least certain and have 
made no attempt to estimate the systemic velocity - the 
source is unresolved and no narrow components are observed 



(the narrowest has FWHM ~ 760 km s~^). In addition to 
the emission line measurements, we have attempted to mea- 
sure the stellar absorption lines (e.g. Ca II K and some of 
the higher order B aimer lines) in a more recently obtained 
VLT spectrum (see lHolt et'allbOO?! ') and estimate a redshift 
of 0.635 ± 0.004. 

Using the above methods, we have confidently determined 
the systemic velocity in 12 of the 14 sources (3C 213.1, 3C 
268.3, 3C 277.1, 4C 32.44, PKS 13454-12, PKS 1549-79, PKS 
1934-63, 3C 459 and with some assumptions, 3C 303.1, PKS 
0023-26, PKS 1306-09 and PKS 0252-71). Due to the issues 
discussed above, the velocities for PKS 1814-63 and PKS 
2135-20 should be used with caution. Further, for all sources, 
the velocity amplitude of the narrow component is coinci- 
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Figure 1. Velocity profiles continued. 



dent with g ravitational motions in an elliptical galaxy (< 
350 km s~ MTadhunter et alJlTosgh . 

The derived/assumed systemic velocities are marked on 
Figure [U and it should be noted that zero on the y-axis of 
the velocity profile plots in Figure [T] does not mark the as- 
sumed systemic velocity, which is denoted by the horizontal 
dashed line. The results are summarised in Table [S] Note, 
whilst the extended gas was used to determine the location 
of the systemic velocity with respect to the various emis- 
sion line components, the redshifts quoted as systemic in 
Table [S] were derived using measurements of several emis- 
sion lines (narrow components) in the nuclear apertures (as 
part of the modelling in Section 4.2) rather than based on 
measurements of [O III] alone. 

In addition to the velocities of the extended emission 
line components in Figure [1] the components of the nuclear 



[O III] lines are shown (see Section 4.2). In the majority of 
the radial velocity profiles, there is a clear offset in veloci- 
ties of the broader components between the extended and 
nuclear apertures of order ~ 100 km s~^. This is most likely 
due to continuum subtraction effects in the nuclear aper- 
tures and/or integrating over a larger aperture. 

HI 21cm absorption is detected in 10 of the 14 sources 
(see Table [6} and the velocities of the detected components 
are plotted on Figure[T] In four sources (3C 213.1, 4C 32.44, 
PKS 1345-1-12 and PKS 1549-79) at least one component 
of HI is in agreement with the assumed systemic velocity. 
However, it is striking that the majority of detected HI com- 
ponents are significantly blueshifted with respect to the op- 
tically derived systemic velocities and, in the case of PKS 
1345+12 at least, appear to trace the optical outflow. In only 
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Figure 2. The effect oo the overall fit to the [O III] AA4959, 5007 
emission line doublet in the nuclear aperture of PKS 0252-71 when 
fitting one (top), two (middle) and three (bottom) Gaussian com- 
ponents to each line in the doublet. In each plot, the overall profile 
(bold line) is overplotted on the extracted spectrum (faint line) 
and the dotted, dashed and dot-dashed lines highlight the rel- 
ative positions and fluxes of the contributing components. It is 
clear that, in the nuclear aperture of PKS 0252-71, a minimum 
of three Gaussian components is required to model the emission 
lines well. 



two sources is the HI redshifted, consistent with infalling gas 
(PKS 0023-26 and PKS 1934-63). 



4.2 Emission line modelling 

As discussed in Section 1, the nuclear emission lines in com- 
pact radio sources are often broad with asymmetric profiles 
requiring multiple Gaussian components to model them. In 
Section 14.11 we identified the rest frame for each source. 
Here, after subtracting the continuum, we model all strong 
emission lines in the nucleus to search for outfiows. In addi- 
tion to multiple Gaussian fitting following H03, we present 
the results of single Gaussian fitting and asymmetry indices. 
These latter two steps will allow us to make comparisons be- 
tween our data and previous work. 

4-2.1 Multiple Gaussian modelling 

The nuclear emission lines in compact radio sources are often 
broad with asymmetric profiles. Hence, single Gaussians are 
unable to adequately model the line profiles and multiple 
components are required. As an example. Figure [2] shows 
how fitting 1, 2 and 3 Gaussians to [O III]AA4959,5007 in 
the nucleus of PKS 0252-71 affects the quality of the overall 
fit. 

When modelling the emission line spectra of AGN, it is 
common to start with the bri ghtest emission lines, e.g. t he [O 
III]AA4959,5007 doublet (e.g. lViUar-Martm et al.lll999t ). The 
[O HI] lines were fitted using three constraints in accordance 
with atomic physics: the fiux ratio between [O III]A4959 and 
[O III]A5 007 was set at 2.9 9:1 (based on the transition prob- 
abilities; lOsterbrockll 19891 ): the widths of the corresponding 
components of each line were forced to be equal; and the 
shifts between the corresponding components of each line 
were fixed to be 48. OA. Note that the fitting program used 
can only work with a wavelength difference and not a ra- 
tio of wavelengths. However, we find the incurred error is 
smaller than our estimated uncertainty. 

Following the technique of H03, we have modelled the 
[O HI] lines using the minimum number of Gaussians re- 
quired to give a good filQ. Hence, for all sources in the sample 
the fits to [O HI] required between 2 and 4 Gaussian com- 
ponents to model the doublet well. The best fitting models 
are presented in Figure[3]and the line data (velocity width^, 
relative shifts) are presented in Tabled The errors quoted in 
Table 4 are a combination of the measurement errors (taken 
from DIPSO) and an estimation of the expected errors based 
on the spectral resolution of the spectra. Typically, for the 
weaker and/or broader components, the dipso error is rep- 
resentative of the true error on the result. For the stronger 
and/or narrower components, the dipso error is very small 
(<<1 per cent) and is not a good indicator of the true er- 
rors. For these sources, we have estimated a percentage error 
based on the spectral resolution and signal-to-noise of the 
data. This is typically of order few- 10 per cent. In addition. 

Note, we define a 'good fit' to the emission lines as that which, 
by eye, has acceptibly low residuals and fits all of the major fea- 
tures (peaks/wings etc.) of the overall line profile well. 
^ Line widths are corrected for instrumental width using the 

equation: FWHMc = ^FWHM^ - IW^ where FWHMc is the 
corrected Full Width at Half Maximum (FWHM) of the line (in 
A), FWHMni is the measured FWHM of the line (in A) and IW 
is the instrumental width or spectral resolution, calculated using 
the measured widths of the night sky lines. 
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the velocities of the nuclear [O III] lines are overplotted on 
the radial velocity profiles in Figure [1] 

From Figure [3] it is clear that a number of sources show 
extreme line widths and shifts, particularly PKS 1345+12, 
4C 32.44, PKS 1549-79 and 3C 459. It is interesting to note 
that these sources account for 4 out of the 5 smallest sources 
in the sample, each having a radio source with linear size <1 
kpc0- For the fifth source in this group, PKS 1814-63, non- 
detection of broader, more blueshifted components may be 
real although difficulties in subtracting the continuum due 
to the presence of a foreground star may also influence the 
result. 

When modelling the emission lines in AGN, it is often 
assumed that one model will reproduce all emission lines. 
This technique has been particularly successful in studies of 
jet-cloud interaction s in pow erful radio galaxies (e.g. Villar- 
Martm et al. 1999; IXavlor e t al. 200^1). Hence, after mod- 
elling [O III] in nuclear aperture of each source, we at- 
tempted to model the other nuclear emission lines with the 
same components and the same velocity widths and shifts as 
the corresponding [O III] model, hereafter the [O III] model. 

As well as the constraints from the [O III] model, some 
emission line doublets required further constraints in ac- 
cordance with atomic physics: the shift between the cor- 
responding components of each line in all doublets (e.g. 
[NeIII]AA3868,3968, [Ne V]AA3346,3425, [N II]AA6548,6583, 
JO I]AA6300,6363 and [S II]AA6716,6731) were set and for 
some doublets (e.g. [Ne III] AA3868,3968, [Ne V]AA3346,3425, 
[N II]AA6548,6583, and [O I]AA6300,6363) the flux ra- 
tios were set based on the transition probabilities. For [S 
II]AA6716,6731, the flux ratio was required to be within the 
range 0.44 < [S II]6716/6731 < 1.42, the ratios correspond- 
ing to the high and low density limits respectively. 

With the exception of PKS 1345-^12 and PKS 1934- 
63, one model (the '[O III] model') fitted the other emission 
lines well. For PKS 1345-1-12, only the narrow component 
was consistently found in all lines whilst in PKS 1934-63, 
none of the emission line components was consistent in all 
lines. However, it should be noted that, due to the complex- 
ity of the emission line profiles in all sources, it was often 
necessary to force the position of the narrow component to 
be at the systemic redshift to obtain a good fit to the lines 
with the [O HI] model. Further, in sources where two narrow 
components are detected, whilst free fitting gave different 
FWHM for the two components, it was also possible to ob- 
tain a good overall fit forcing the two narrow components to 
have equal FWHM (e.g. 3C 268.3, 4C 32.44, PKS 0252-71, 
PKS 1306-09, as discussed in Section [4. l|l . In these sources 
the derived systemic redshift was an average of the redshifts 
of the two components. 

In addition to emission originating from the NLR, we 
have detected broad components to the permitted lines (due 
to the signal-to-noise of the data, often only detected in Ha), 
i.e. from the BLR, in three sources: 3C 268.3, 3C 277.1 and 
PKS 1549-79 (H06). 3C 277.1 is known to be a quasar and 
so the detection of a broad component to the permitted 
lines was not surprising. However, the detection of a broad 



° Note, for 3C 459, whilst the overall radio source is highly ex- 
tended the compact core is consistent with the other compact 
radio sources in the sample. 



component in Ha in 3C 268.3 has led to us re-classifying 
this radio source as a BLRG rather than an NLRG. The 
detection of broad Ha and Paa in PKS 1549-79 was reported 
in H06. The line data for the BLR components is presented 
in Table H 



4.3 Single Gaussian Modelling 

Whilst multiple Gaussian modelling has worked well for the 
majority of sources in the sample, two GPS sources (PKS 
1345+12 and PKS 1934-46) were notable exceptions. In 
PKS 1345+12, a common narrow component was observed 
in all lines although different emission lines required differ- 
ent velocity widths and shifts for the broader components 
to provide good fits (H03). In PKS 1934-46, the situation 
is somewhat worse - it is less clear that a common nar- 
row component exists for all lines when observed at this 
resolution / signal-to-noise. 

For PKS 1345+12, H03 suggested a stratified ISM with 
three distinct regions, each responsible for the emission of 
one of the kinematic components. Hence, the different inter- 
mediate and broad components could be explained by gradi- 
ents in ionisation potential and/or critical density across the 
regions emitting the intermediate and broad components. 
Stratification of the NLR was suggested as the reason for 
correlations between emission line FWHM (when fitting a 
single Gaussian) an d ionisation potential or cr itical density 
in Seyfert galaxies (Ide Robertis fc Osterbrock| [l984. 198d). 
This is consistent with models in which there is a continuous 
variation of density, ionisation and velocity across a spatially 
unresolved inner narrow line region (INLR). This te chnique 
was applied to a sample of extended radio galaxies bv lTavlon 
(2004). In the latter work, no significant correlations were 
found in Cygnus A or in the sample of NLRGs, but 3/4 
BLRGs studied did show significant correlations with ioni- 
sation potential and/or critical density, a feature expected 
for an INLR on a scale less than that of the central obscuring 
torus (r < 100 pc). 

We have modelled all nuclear emission lines in all 
sources using single Gaussian free-fits. No clear correlations 
were observed for the majority of th e sourc e s, in cluding the 
BLRGs, in contrast to the results of'Taylor" ("2004). Figure H 
shows plots of the single Gaussian FWHM versus both criti- 
cal density and ionisation potential for the smaller sources in 
the sample: PKS 1345+12, PKS 1549-79 and PKS 1934-63. 
Whilst visual inspection of the plots suggest trends, par- 
ticularly in critical density. Spearman Rank analysis shows 
these visual trends to be insignificant. 

The line widths (FWHM) derived from the single Gaus- 
sian modelling can also be used to compare this sample of 
compact radio sources to other samples of radio sources 
presented in the literature. The measured single Gaussian 
FWHM of the nuclear [O HI] lines for this sample are pre- 
sented in Tableland are compared to a sample of extended 
radio sources in Section 5. 
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Figure 3. Models to the [O III] AA49.59,. 5007 emission line doublet in the nuclear apertures of all 14 sources in the sample. The faint solid 
line represents the continuum subtracted spectrum and the bold solid line is the overall model fit to the doublet. The components are: 
narrow (dotted), intermediate (dashed), broad (dot-dashed) and very broad (dot-dot-dashed) for all sources except 3C 277.1 where the 
dot-dot-dashed line traces the BLR component of H/3 (see Section 3.3 for the kinematic definitions). All plots are on the same wavelength 
(i.e. velocity) scale to allow for comparisons between the outflows. 



5 DISCUSSION 

5.1 Extreme emission line outflows 

The nuclear emission lines in all of the compact radio sources 
in this sample are liigiily broadened, with complex emis- 
sion line profiles, and require multiple Gaussian components 
to model them. These components are often shifted signifi- 
cantly with respect to one another and are therefore likely 
to trace flows in the nuclear emission line gas. However, in 



order to distinguish between inflows and outflows, it is nec- 
essary to accurately define the systemic velocity. 

In Section 4.1, following a similar approach to that used 
for PKS 1345-hl2 (H03) and PKS 1549-79 (H06), we argued 
that the extended narrow component (s) represents the am- 
bient, quiescent ISM in the galaxy halo. We are confident of 
this result in 12 of the 14 sources studied (see Section 4.1) - 
kinematically the narrow component is the least disturbed 
component and is typically the only extended component. 
In some of the sources in which HI absorption is detected. 
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Figure 3. continued. Models to the [O III] AA4959, 5007 emission line doublet. 
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Table 4. Rest frame emission line modelling results. Columns are: (a) object name; (fe) component following the definitions above: n 
(narrow; rn and bn denote the 'reddest' and 'bluest' component where two narrow components are detected), i (intermediate), b (broad), 
vb (very broad) and BLR (broad line region component; i.e. a highly broadened component only observed in the permitted lines); (c) & 
(d) velocity width (FWHM) and error of each component; (e) & (/) velocity shift and error of the component from the systemic velocity. 
Note, the results for PKS 1345+12 and PKS 1934-63 are for [O III]. In sources for which free fitting gave two narrow components with 
different widths but the lines can also be modelled adequately by forcing the FWHM of the two narrow components to be equal, the 
FWHM of this alternative model for the narrow component is given and marked f. 



Object Velocity A Velocity A 
width shift 

km s~^ km s"-"^ km s~^ km s" 

(a) (b) (c) (d) (e) (/) 



3C 213.1 
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the velocity of the narrow component is also consistent with 
a deep, narrow component of HI w hich could be associat ed 
with a circumnuclear disk or torus (jMorganti et aLlbOOll ). 

Previous studies have often struggled to determine the 
relationship between the HI, the optical kinematics and the 
systemic velocity (i.e. is the HI at rest, or tracing a flow 
of gas?) due t o the inaccuracies in tying down the systemic 
velocities (see lMorganti et al.ll200l] for a discussion). In this 
paper, we have accurately established the systemic velocities 
in 12 out of 14 sources using detailed optical emission line 



modelling and have compared them to the HI data available 
in the literature in Figure [1] and Table |6] 

HI absorption is detected in 10/14 sources. Compar- 
ing the optical and HI data, it is striking that the major- 
ity of absorption components are significantly blueshifted 
with respect to the optical systemic velocity, often by sev- 
eral hundreds of km s~^, and represent outflows of neutral 
gas. Only three sources have HI components consistent with 
the optical systemic velocity (3C 213.1, 3C 268.3, and PKS 
1345+12) which may be consistent with a circumnuclear disk 
or torus. However, in PKS 1345+12, the latest data show the 
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Table 5. Various emission line pajrameters used in the statistical analysis. Columns are: (a) object; (6) heliocentric redshift of the systemic 
velocity; (c) maximum shift velocity taken as the shift between the broadest NLR component and the systemic velocity (km s~^); d 
FWHM of the [O III]A5007 line when fitting a single Gaussian (km s~^); (e) the asymmetry index, AI20 and (/) the kurtosis parameter, 
R20/50- * indicates total range of HI absorption and the PKS 1306-09 data is from Raffaella Morganti (priv. comm.). f assumed optical 
z. " The shift between the two detected narrow components. ★ Redshift derived from the stellar absorption lines. 
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Figure 4. Correlation plots for the nuclear apertures of PKS 1345+12, PKS 1549-79 and PKS 1934-63. For each source the plots arc: 
top : rest frame FWHM (km s-^) versus log(critical density cm ^) and bottom : rest frame FWHM (km s ^) versus lonisation potential 
(eV). 
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component of HI at the systemic velocity i s actually most 
likely associated with an off-nucleus cloud (jMorganti et alj 
I2OO4I ). Unlike the optical outflows, both the narrower and 
broader line components have signiflcant blueshifts. Multiple 
HI components are observed in 5 sources. In PKS 1345+12, 
these are broadly consistent (velocity width and shift) with 
the optical components. Finally, two sources contain com- 
ponents consistent with infalling clouds (PKS 1934-63 and 
PKS 0023-26). 

One of the key trends highlighted above (Section 4.2) 
is the detection of significantly shifted emission line com- 
ponents observed in the majority of sources in the sample. 
These are predominantly blueshtfted, with the broader com- 
ponents blueshifted by the largest amounts. Figure [7] shows 
the distribution of shiftfl for the sample. All sources except 
PKS 1306-09 - which has no 'broad' component, only two 
narrow components - are included in this plot and the shifts 
used are for the NLR components only. 

For all but two sources (PKS 0023-26 and PKS 0252- 
71) the broad component is blueshifted with respect to the 
systemic velocity with velocities up to 2000 km s~^ although 
most sources occupy the range -900 < Ushift < 100 km s~^. 
Figure [7] further sub-divides the sample into the larger CSS 
sources and the smaller GPS sources. The effects of radio 
source size and orientation are discussed below. 

The [O III] line emission is known to be resolved in 
three of the sources in this sample: 3C 268.3, 3C 277.1 
and 3C 303.1. Thes e sources have been imaged at high res- 
olutio n w ith HST llde Vries et al] 1 19971 . Il999l : lAxon et al.1 
I2OOOI ) aud io 'Pea et al.l (|2002l ) have obtained corresponding 
HST/STIS spectroscopy of 3C 277.1 and 3C 303.1. For these 
sources, virtually all of t he line emission de tected by both 
Ide Vries etHI (|l999h and I Axon et al.l (|2000l ) hes within our 
slit. In our 2-D spectra, we observe clear extended line emis- 
sion in 3C 268.3 and 3C 277.1 although the line emission im- 
aged by HST lies within our defined nuclear aperture, within 
which no obvious spatial structure is resolved. The extended 
line emission that we detect will be presented in Holt et al. 
(2008, in prep.) along with a discussion of extended aper- 
tures. It is difficult for us to c ompare our kinematic results 
to those of lO'Dea etHI (I2OO2I ) as whilst their spectra are of 
significantly higher spatial and spectral resolution compared 
to ours, from the figures in their paper, the lines are at much 
lower signal to noise than our spectra and they have, for the 
main part, modelled the emission lines with single Gaus- 
sians. Depending on the dominance of different components, 
this can lead to erroneous line shifts (see e.g. HOG). Com- 
paring our single Gaussian modelling, we observe similar line 
widths (FWHM~500km s"^) in 3C 277.1 but significantly 
higher widths in 3C 303.1 (F^yHM^ 800km s"^ compared 
to ~500km s"^ in lO'Dea etal]l2002h . Further, the resolved 
[O HI] emission may account for the narrow line splitting in 
3C 268.3. 

From the results presented here, it is clear that, with the 



^ Note, for the sample described in this paper, the broadest (i.e. 
broadest NLR component) to systemic velocity shift is used rather 
than the broadest to narrow velocity shift. This is because sev- 
eral sources have multiple narrow components which are believed 
to trace rotation curves in the galaxy rest frame and using the 
systemic velocity is the most consistent method for all sources in 
the sample. 



exception of PKS 0023-26, PKS 0252-71 and PKS 1306-09, 
all sources in the sample have evidence for fast, blueshifted 
fiows in the circumnuclear ISM. However, from the kinemat- 
ical evidence alone it is impossible to distinguish between 
material in outflow on the side of the nucleus closest to the 
observer and mat erial infalling on the far size of the nucleus. 
iHolt et al.l (|2003l ') argued that, for PKS 1345+12, as the red- 
dening increased significantly with FWHM, the blueshifted 
components were likely to be consistent with material being 
observed on the observer's side of the nucleus and therefore 
tracing an outflow. We will discuss the issue of reddening in 
a future paper (Holt et al. 2008, in prep.). 



5.2 The role of orientation 

The original wor k on the compact fiat-sp ectrum radio source 
PKS 1549-79 bv'Tadhu nter et al.l l|200ll ) suggested that the 
direction of jet propagation was oriented close to the ob- 
server's line of sight (LOS) as: 

(i) PKS 1549-79 has a flat radio spectrum, often associ- 
ated with radio loud quasars whose axes are close to the 
observer's LOS; 

(ii) PKS 1549-79 has core-jet radio morphology. As radio 
galaxies have two-sided radio jets, observing one-sided ra- 
dio structure is usually interpreted as an orientation effect 
with the observer's LOS aligned close to the direction of jet 
propagation (see H06, Figure lojfl- 

Simil arly, radio maps of PKS 1345+12 (e.g. iLister et al.l 
I2OO3') show highly asymmetric radio jets, again suggesting 
the radio source is pointing close to the observer's LOS. This 
argument is supported by the det ection of a broad (FW HM 
~2600 km s~^) component Paa jVeilleux et al.lll997h and 
the detecton of a point source comp onent in high-resolution 
near-IR images (|e vans et al]|l999|i. Howev er, this orienta- 



tion is disputed bv lLister et all l|2002l. |2003 ) - see H03 for a 
discussion. 

Interestingly, PKS 1345+12 and PKS 1549-79 contain 
two of the most extreme outflows in the sample discussed in 
this paper (~2000 km s~^ and ~680 km s~^ respectively). 
Perhaps, then, the particular orientation of these objects 
with respect to our LOS is the reason why we observe such 
extreme kinematics in these sources. In contrast, for exam- 
ple, PKS 1934-63 has a weak core and is a symmetric double, 
as well as being one of the most compact sources in the sam- 
ple, and has one of the smallest observed outflow velocities. 

Due to the angular scales of the compact radio 
sources, compact flat-spectrum radio cores for the sources 
in thi s sample are not gene rally detected, even with VLBI 
(e.g. iTzioumis et"aLl I2OO2I ) - the core is only detected 
in 5/14 sources. It is therefore not possible to deter- 
mine the radio source orientation accurately for most of 
the sample using the standard meth ods, for example, the 
R (core-dominance) p arameter (e.g. lOrr fc Brownd Il982l : 
ICiovannini et al]l200ll l , radio flux density comp arison of the 
jet and counter-jet (e.g. ICiovannini et al.|[200ll ). jet proper 



° It should be noted that highly asymmetric radio structures are 
also consistent with jet-cloud interactions in a non-homogeneous 
ISM. 
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Table 6. Summary of the HI data for the sample compiled from the literature. Columns are: 1) object, 2) & 3) HI velocity shift and error 
with respect to the optical systemic velocity derived in this paper where negative velocities imply blueshifted HI; 4) & 5) HI full width at 
half maximum (FWHM) and estimated uncertainty (° full width at zero intensity (FWZI) of the broad, shallow absorptio n); 6) optical 
depth qu oted fractional absorpt ion); 7 ) HI column density d erive d assuming Tgj,i„ = 100 K an d 8) reference s : vO: Vcr on-Cettv et al. 
1 20001 ): m l: iMorganti et al.l feoOll ); m3: iMorganti et al.l (|2003| ): v3: IVermeulen et al.l ((200.3); g6: iGupta et all I2OO6I) : H06: iHolt et al. 
1 20061): 16:lLabi ano et aP l|200d ): m: Morganti, private communication. * component of HI consistent with the optical systemic velocity, 
t IVeron-Cettv et al. 1 1I2OOOI) " give e.g. FWZI and we have estimated the FWHM from their figures but also give their data here. 
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motion measurements ( e.g. Giovannini et al] l200lh or jet 



motion modelling (e.g. iLister et al. 20031 ). Three sources 
have quoted orientation inf ormation in the literat ure: PKS 
1549-79 {R2.3GH. 



(Ri.&GHz = 5.65; _ 
IMorganti et al.|[200 



1.310: IMorganti et al.ll200ll ). 3C 459 
Morganti et all 1 19931 : ?.nnr = 0.394; 



11) 

elling by Lister et al.ll200 j : 82° to LOS with opening angle of 
46°, although this result is inconsistent with the arguments 
presented in H03). 

Despite the lack of radio information to accurately de- 
termine the radio source orientations, it is still interesting to 
make a rough estimate of the importance of orientation with 
respect to the observed outflow velocities. For the remainder 
of the sample, we have therefore adopted the rather crude 
approach of comparing the relative extents of the radio jets 
on either side of the nuclei in the radio maps. Assuming 
both jets are generated by the same central source, if they 
expand through empty space they should have similar in- 
trinsic extents on either side of the core and any observed 
differences will be due to orientation effects. It should be 
noted that, should the jets expand through a dense ISM, 
interactions between the radio source and this ISM could 
significantly alter the path of at least one of the jets. Hence, 
a highly asymmetric ISM can alter the relative extents of 
the jets and mimic the effects of orientation with respect 
to the observer and these results should therefore be used 



and PKS 1345-1-12 (jet motion mod- 



with caution. However, iBest et all (Il995l ) studied the angu- 
lar asymmetries in a sample of extended FR II 3CR sources 
and found the distribution of asymmetry angles to be con- 
sistent with the predictions of unified schemes. 

We have therefore classified the sources into three broad 
categories: 

(i) Close to the LOS. All sources with an obvious core-jet 
and a large core/extended radio flux ratio e.g. PKS 1549-79, 
PKS 1345-H2 and 3C 303.1 

(ii) Close to the plane of the sky. Symmetric radio mor- 
phology (e.g. similar jet extents and/or fluxes): e.g. PKS 
1934-62, PKS 0252-71 and PKS 0023-26. 

(iii) Intermediate. Sources not consistent with the other 
two categories. 

Figure [5] shows the grouping with respect to the largest out- 
flow velocity. 

Despite the crudeness of the method, the data suggest 
that the orientation of the radio source to the observer's 
line of sight/the asymmetry of the radio morphology may 
be an important factor in determining the maximum out- 
flow velocity observed. Sources with jets most likely point- 
ing towards us contain the most extreme outflow velocities, 
and sources likely to be close the plane of the sky exhibit 
the least evidence for fast outflows in the NLR. Only two 
sources in the sample appear to be 'mis-fits' in this classifi- 
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Table 7. Comparison of the maximum outflow velocity with the 
largest projected linear size of the sources in this sample. " com- 
pact core. 



D < 1 kpc D > 1 kpc 



quasars. Whilst the BAL outflows are on a different scale 
to those observed in compact radio sources (on the scale of 
the BLR rather than the NLR), it has been suggested that 
the extreme velocities observed may be partly due to an 
orientation effect (e.g. IWevmann et al-lfigQll : lElvij|2000l ). 
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Figure 5. Histogram showing the putative orientation of the 
sources with respect to their outflowing velocities. 



cation scheme - 3C 277.1 and 4C 32.44. This is most likely 
due to both the over-simplification of the method and inher- 
ent scatter in the data. Note also, whilst the flat-spectrum 
core appears strong in the radio map of 3C 277.1 and we ob- 
serve strong optical signatures of the quasar in our optical 
spectra (broad permitted lines and a strong blue power-law 
continuum), the radio emission is resolved into two distinct 
lobes which suggest the radio axis is not closely aligned with 
our line of sight. We have therefore placed 3C 277.1 in the 
'intermediate' category. 

Such a strong trend in orientation-outflow velocity has 
been suggested for other classes of AGN. The most extreme 
outflows observed in any type of AGN are those in BAL 



5.3 Outflow driving mechanism 

It is clear that extreme emission line kinematics exist in ah 
of the compact radio sources in this sample. The emission 
line flows are likely to trace outflows in the emission line gas, 
similar to the results presented for PKS 1345+12 (reddening; 
H03) and PKS 1549-79 (source orientation; H06). 

If compact radio sources are truly young AGN in which 
the outflows have not yet swept aside the natal cocoon, then 
because of the large gas and dust concentrations in the cen- 
tral regions, the effects of all types of outflow (quasar winds, 
jets, starburst superwinds) are likely to be more visible as 
there is more warm/cool gas around for the winds and/or 
jets to interact with. Whilst jet-cloud interactions are ex- 
pected to produce extreme kinematics, all of these scenarios 
can potentially explain broad line widths and large velocity 
shifts and so, from the kinematical evidence alone, it is not 
possible to distinguish between the different outflow driving 
mechanisms. 

Some of the best evidence for jet-driven outflows is pro- 
vided by the co-alignment and similar scales of the radio 
source and the optical emission line structures. More than 
30 CSS sources have been observed with HST in both broad 
an d narrow bands (to isolate [O II] and/or [O I I I] emi ssion) 
bv lde Vries et all (|l997. , 1999) and lAxon etal] l|200d ). The 
de Vries et al. samples include 3C 213.1, 3C 268.3, 3C 277.1 
and 3C 303.1 and the Axon et al. sample includes 3C 268.3, 
3C 277.1, 30 303.1 and PKS 1345+12. These imaging stud- 
ies reveal that the optical and radio emission is both: 

(i) on similar scales (30-90% of the optical line emission 
is concentrated within a few kpc and the radio sources are 
either completely embedded within this optical line emitting 
gas or extend just beyond it; lAxon et al.ll2000l sample) and; 

(ii) strongly aligned (typically < 10° in all sources, in- 
cluding the sources also found in our sample, in which the 
data could be accurately registered in the de Vries et al. 
sample, and in 6/11 sources in the Axon et al. sample, often 
elongated into jet-like structures). 

The radio-optical alignment is also observed across the en- 
tire redshift range probed by these samples (0.1 < z < 
1.5), rather than conflned to highe r redshifts as in more 
extended radio sources (z > 0.6; e.g. iMcCarthv et al.|[l987l : 
Ide Koff et al.|[l996l '). Such close association between the opti- 
cal emission line gas and the radio source suggest that the ra- 
dio source is strongly interactin g with the amb ient medium 
as it expands through it ifO'Dea et al.l |2002| ). Follow-up 
HST/STIS spectroscopy of three CSS sources with resolved 
[O III] line emission (including 3C 277.1 and 3C 303.1) pro- 
vides further evidence that the outflows are likely to be 
driven by the expanding radio jets; the emission lines have 
complex, broad proflles (FWHM ~ 500km s~^) which are 
offset with respect to the systemic velocity b y 300-500km 
s~^ in the region(s) of the radio hotspots l|0'Dea et al.l 
and the line ratios are consistent with a mixture of 
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fast shocks (500-lOOO km s ^) and photoionisation/precursor 
l|Labiano et al.ll2005h . 

Finally, PKS 1345+12 and PKS 1549-79 were re- 
cently observed at h igher resolution with HST/ACS 
l|Batcheldor et al.ll2007l ). In these sources, the optical line 
emission is clearly concentrated in the central regions, ruling 
out galaxy wide starbursts as the dominant outflow driving 
mechanism. Further, the radio emission in both sources is 
on similar scales to the optical line emission. However, the 
region of line emitting gas is only marginally resolved and, 
whilst there is a suggestion that the isophotes may be elon- 
gated in the direction of the radio axis, this evidence is far 
from conclusive. 

Despite the strong evidence for radio/optical alignments 
and the co-spatial scales of the optical and radio emission, we 
find no evidence for anti-correlations between the emission 
line widths and the ionisation potential (Figure that are 
a clear signature of shocks in some extended radio sources. 
Hence, to confidently distinguish between the different ioni- 
sation mechanisms, it is necessary to combine the kinemati- 
cal and imaging data with a detailed study of the line ratios 
(Holt et al., 2008, in prep.). 



e 10 



12 — 

S 8 - 
o 

067 
I 4 - 

2 r 
^ 



Extended Radio sources (D > 30 kpc) 



1000 1500 2000 

FWIIM (krn/s) 



Compact Radio sources (D < 30 kpc) 



1000 1500 

J'ft'HM (km/s) 



5.4 Are the kinematics in compact radio sources 
more extreme than in extended radio sources? 

As discussed in Section 1, previous data have suggested 
that the kinematical properties of the nuclear emission 
line gas are different in compact radio sources compared 
to more extended radio sources. For example, comparison 
of the emission line profile s in samples of co mpact (e.g. 
Gelderman k. Whittlelll994l ') and extended (e.g. iBrothertonI 



19961 ) radio sources in the literature suggest that the line 



profiles of the compact radio sources are often broader and 
more asymmetric. However, to date, no attempt has been 
made to properly quantify the differences between the emis- 
sion line profiles of compact and extended radio sources. 

Many of the spectra in the literature are not of sufficient 
signal-to- noise and/or resolution to model the emission line 
profiles in detail and so we have used three line profile pa- 
rameters on the strong [O HIJA5007 line to compare our 
sample to the data in the literature, for example, i) the line 
width (FWHM) derived from single Gaussian modelling ; 
ii) the asymmetry index, AI20, and iii) the kurtosis param- 
eter, R 20/50. The latter two are defined in iHeckman et aP 
l|l98ll ). All the results for our sample are presented in Ta- 
ble [5] As our sample of compact radio sources is small, we 
have formed a larger sample including data from the lit- 
erature (see Table IH]). Our comparison sample of extended 
radio sources for these three pa rameters compri s es all radio 
loud object s in th e samples of iHeckman et ahl (jl984l l and 
IBrothertonI (j 19961 ) (see Table [8|. The latt er sample com- 

J rises solely quasars whereas the sample of i Heckman et all 
19841 ) includes both radio galaxies and quasars. 

Figure [6] shows that the distribution of FWHM of [O 
III]A5007 is markedly different for the compact and ex- 
tended radio source samples. For the extended sources, 
the distribution peaks at ~ 200-300 km s~^ which is con- 
sistent with gravitatio nal motions in an elliptical galaxy 
ijTadhunter et al.ll 19891 ). However, more extreme kinematics 
are observed in the sample of compact radio sources - the 
spread is larger and it peaks at FWHM ~1000 km s~^ - ve- 



Figure 6. Histograms showing the distribution in observed op- 
tical emission line FWHM of [O III] in sa mples of extended 
(D > 30 kpc; from IHeckman et al.|[l984l and IBrothertonI Il996l) 
and compact (D < 30 kp c; this sample plus sources from 
Gelderman fc Whittle! lll994l) . four sources from the sample of 



Hirst et al.l ||2003|) and PKS 1718-49 from Fosb urv et al.l l|1977:)) 
The shaded part of the bottom plot highlights the sources from 
this sample. 
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Figure 7. Histograms showing the maximum nuclear outflow 
velocity for this sample of compact radio sources (bottom) 
and a c ompa rison sample of extended radio sources taken from 
iTavloJ (Eool). Note, for the extended sources the shift plotted 
is broadest— narrowest component whilst for the compact sources, 
the shift between the broadest-systemic velocity is used as in 
a number of sources, two narrow components are observed (see 
Section 4.1). For the compact radio sources, the plot is further 
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Table 8. Parameters for the comparison samples used in both this paper and in Paper II. Columns are: (a) sample name; (b) redshift 
range; (c) radio power range at 5 GHz except " at 178 MHz, ^ at 1.4 GHz and no radio data was given for this sample; (d) number of 
sources used; (e) completeness of the sample, f see Section 2 for details. Note, for sources which are included in more than one sample, 
data are included only once and are taken from the sample with the highest quality data. * This comparison sample is only used in Holt 
et al., 2008, in prep. 



Sample 



(a) 



Radio power N completeness 
log PsGHz 
(W Hz-i) 

(c) id) (e) 



Compact radio sources 

This paper 

Gelderman fc Whittle (1994) 
Hirst et al. (2003') 
PKS 1151-34 
PKS 1718-643 



< 0.7 

< 0.9 
0.5-3.6 

0.258 
0.014 



26-28 
22-29 
> 26.5" 
26.8 
24.3 



14 
16 
9 



completef 
not complete 
complete 



Extended radio sources 
2 Jy 

radio loud objects from jIeckmar^_et_aL Q_984 
radio loud objects from ^rotherton (_1996) 
Taylor (2001) 




25- 29 
> 24.5'' 

C 

26- 29 



36 
45 
60 
12 



complete * 
not complete 
not complete 
complete 



locities which are too large to be explained by gravitational 
motions unless the gas was located close to the central black 
hole (i.e. in the BLR or an INLR). A one-tailed Kolmogorov- 
Smirnov (K-S) test shows the distributions are different at 
the 99.9 per cent confidence levelf]. Whilst this effect might 
be due, in part, to orientation effects in both compact radio 
sources (see Figure 5) and/or foreshortening of larger radio 
sources, it has been argued that the majority of CSS sources 
can not be larger radio s ources foreshortened by projection 
effects (|Fanti et aLlliggd ). 

Whilst the line widths show a clear difference between 
the compact and extended sources, a K-S test shows that 
neither the asymmetry index or the kurtosis parameters are 
significantly different. Whilst we may initially expect the 
broad, blueshifted components to give highly asymmetric 
line profiles, examination of Figure [3] shows in many sources, 
either the peak fiux of the broadest components is often 
below the 20% fiux level (i.e. asymmetries would only be 
evident in parameters using e.g. the 10% fiux level as in PKS 
1345+12), or the emission lines appear symmetric although 
the dominant component may be highly shifted (e.g. PKS 
1549-79). 

In addition to the line widths/profiles comparison to 
older data in the literature, we have compared our sample 
to a recent study of a sample of extended 3C radio sou rces 
taken with similar observational set-up l|Tavloij |2004| ) en- 
abling comparisons between, for example, the shifts between 
the different components of a line (see Table [S]). A histogram 
of the broadest-narrow velocity shifts for both samples in 
shown in Figure [T] The distributions are markedly different 
(99.9%: K-S test) with the compact radio sources spanning 
a much larger range in velocity (out to '^2000km s~^ com- 



^ Whilst the K-S test reveals a strong trend, the reader should 
also remember that the various samples have different selection 
criteria which may affect the test to some degree. 



pared to ~600km s"'^). It is interesting to note that this 
trend of more extreme velocities in smaller sources is also 
evident within the sample of compact radio sources, with 
the two most extreme velocities being detected in smaller 
GPS sources (PKS 1345-hl2 and 4C 32.44). Indeed, aU but 
one of the most extreme velocities (v > 500 km s~^) are de- 
tected in sources with projected linear size D < Ikpc whilst 
the cluster of sources in the low outflow velocity region (-100 
< Vshift < 100 km s~^) are predominantly the larger CSS 
sources (Figure[7]and TableO. Such velocity evolution with 
source size is consistent with the idea that compact radio 
sources are young, starting as GPS sources and expanding 
to become CSS then extended radio sources. 

Whilst the nuclear regions of extended radio sources 
appear to be less extreme than in more compact ra- 
dio sources, extreme velocit ies (v > 1000 km s~^; e.g. 
ISolorzano-Inarrea et al-lboOll ) are commonly observed in ex- 
tended radio sources, but in the extended emission line re- 
gions (EELRs). EELRs are typically on similar scales as, 
and are aligned with, the radio axis and are therefore asso- 
ciated with jet-cloud interactions. Hence, in extended radio 
sources, the regions emitting the broader components are 
often significantly displaced from the nuclear aperture, co- 
incident with the radio emission. 

With this in mind, it is not surprising that extreme 
kinematics are observed in the nuclear apertures of com- 
pact radio sources. Compact radio sources are on the scale 
of the host galaxy (<15 kpc), often contained within the 
nuclear regions (<1 kpc), and so the signatures of inter- 
actions between the radio jets and the ISM will occur in 
the circumnuclear regions. Indeed, it is in the nuclear aper- 
tures of compact radio sources that we observe both broad, 
blueshifted components, which are generally spatially unre- 
solved, as well as any quiescent components similar to those 
observed in the nuclear regions of extended radio sourcces. 
The consistency between the spatial scales of the emission 
line regions and the radio emission has been confirmed in 
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six sources in our sample (3C 213.1, 3C 268.3, 3C 277.1, 3C 
303.1, PKS 1345+12 and PKS 1549-79) using HST broad- 
and narrow-band imaging Vries et al ] |l997l : Ia xon et al.l 
[2000 ; Bat chcldor ct al. 200^ Hence, whilst the kinematics 
observed in the nuclear regions of compact radio sources 
are significantly different to those observed in the nuclei of 
more extended radio sources, they are entirely consistent 
with compact radio sources being a young, small scale ver- 
sion of the jet-cloud interactions observed n the extended 
sources with aligned radio and optical emission. 



6 SUMMARY, CONCLUSIONS AND FUTURE 
WORK 

It is clear that all compact radio sources show evidence for 
disturbed kinematics, with large line widths and shifts with 
respect to the galaxy rest frame in the optical emission lines. 
The main conclusions of this paper are: 

• The extended emission line halo. For the majority 
of sources (12/14), we were able to determine the systemic 
velocity with confidence. For consistency, all techniques fo- 
cused on the extended narrow component emission. In two 
sources, a smooth 'rotation curve' was observed, but in most 
sources, the extended narrow component emission was either 
not settled, as in PKS 1345-1-12, or split into two narrow 
components, assumed to represent a rotation curve in which 
the central regions were unresolved although this may also 
be a signature of bi-polar outflows. 

• Extreme emission line outflows. All but three 
sources show evidence for outflows in the circumnuclear 
ISM. The most extreme outflow is in the GPS source PKS 
1345+12 (~2000 km s~^). Interestingly, the second most 
extreme outflow was observed in another GPS source, 4C 
32.44. As well as radio source size (CSS or GPS), the orien- 
tation of the radio source to the observer's line of sight may 
also be important, with higher outflow velocities observed 
in sources pointing towards the observer. 

• Blueshifted HI. HI absorption is detected in 10/14 
sources, with multiple components observed in 5 sources. 
The majority of HI components (narrow and broad) are sig- 
nificantly blueshifted with respect to the systemic velocity 
and trace outflows in the neutral gas. In PKS 1345+12, (and 
others?) the outflowing HI components are broadly consis- 
tent with the emission line components. In only two sources 
(PKS 0023-26 and PKS 1934-63) is the HI redshifted, con- 
sistent with infalling gas. 

• Kinematical evidence for shocks. As well as the ex- 
treme line splitting/outflows, the emission line components 
are also highly broadened. Again, the two most extreme 
FWHM in the NLR gas are observed in GPS sources: PKS 
1345+12 (~2000 km s~^) and 4C 32.44 (~3500 km s~^). 
Highly broadened components are almost exclusively con- 
fined to the nuclear apertures and therefore on the scale of 
the radio source. Higher resolution HST studies have re- 
solved the [O HI] emission line gas in several sources in 
this sample, revealing that the emission line gas is both on 
the same scale as, and is strongly aligned with, the radio 
source. This is consistent with observations of high redshift 
extended radio sources in which the broadest components 
are observed in regions coincident with the radio emission 



(e.g. ISolorzano-Ifiarrea et"al] I2OO1I ). The suggested impor- 
tance of orientation on the observed outfiow velocity also 
suggests the acceleration is confined to a spatially small re- 
gion (i.e. coincident with radio jets) rather than ocurring 
across the entire nuclear region (i.e. due to a quasar- or 
starburst wind). 

• Smaller source — more extreme kinematics. These 
results provide further evidence that the radio source size 
is also important in determining the outflow velocity and 
the dominance of shocks. Shocks, and the kinematics as- 
sociated with them, are predom inantly con flncd to smaller 
extended sources (D < 120 kpc: iBest et al.l 2000). The sta- 
tistical tests presented in this paper show that the nuclear 
kinematics, in particular the emission line shifts, are more 
extreme in compact radio sources than in the nuclear regions 
of their extended counterparts at the 99.9% confidence level. 
Similarly, the most extreme outfiows in the compact radio 
sources occur in sources with the smallest projected linear 
sizes, generally with D < 1 kpc, the GPS sources. Whilst this 
effect might be due, in part, to foreshortening, it has been 
argued that the majority of CSS sources can not be la rger ra- 
dio s ources foreshortened by projection effects ()Fanti et al.l 
Il990h . 
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